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Abstract. We explore the dependence of the amplitude of stellar dynamo cycle variability (as seen in the Mount Wilson Ca 
II HK timeseries data) on other stellar parameters. We find that the fractional cycle amplitude A cyc (i.e. the ratio of the peak- 
to-peak; variation to the average) decreases somewhat with mean activity, increases with decreasing effective temperature, 



but is not correlated with inverse Rossby number Ro 
frequencies u) cyc /£l along two, nearly parallel branches. 



We find that A cyc increases with the ratio of cycle and rotational 
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1. Introduction 

In a recent series of papers (Brandenburg, Saar, & Turpin 
1998; Saar & Brandenburg 1999 [=SB], 2001), we have been 
exploring relationships between magnetic cycle periods P cy c, 
rotation periods P ro t, and other stellar properties. By combin- 
ing P cyc obtained from observations of Ca II emission (e.g., 
Baliunas et al. 1995 [=Bea95]), photometry, and P ro t vari- 
ation (e.g., Lanza & Rodono 1999), we found evidence for 
trends between cycle (w C yc) and rotational (il) frequencies 
(Saar & Brandenburg 2001), and between u> cyc /il and both 
i?^ K and Ro' 1 (Brandenburg et al. 1998; SB). Here, Ro- 1 
is the inverse Rossby number and R' HK is the Ca II HK flux, 
corrected for photospheric contributions and normalized by 
the bolometric flux (see Noyes et al. 1984). 

Another important, but less studied observational prop- 
erty of stellar dynamos is the cycle amplitudes (i.e., the mean 
amplitude of the cyclic variability as seen in some activity 
diagnostic). There is a relationship between cycle amplitude 
and period seen in the Sun, where shorter cycles tend to be 
followed by stronger ones (e.g., Hathaway, Wilson & Reich- 
mann 1994), but there has been little work along these lines 
in stars. Focusing on inactive stars, Soon et al. (1994) sug- 
gested that the fractional peak-to-peak cycle amplitude seen 
in chromospheric Ca II HK emission, 

A cyc = AR' HK / (i?H K ), 
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decreased linearly with \og(£l / uj cyc ) 2 . Studying a larger sam- 
ple, Baliunas et al. (1996) found that among inactive stars, 
A cyc oc (w C y C /f2) - 9 . Results for active stars were less clear. 

Baliunas et al. (1996) analyzed all cycles detected by 
Bea95, independent of their quality. In our earlier work on 
cycle periods we found it very useful to begin by taking only 
a critically chosen sample of reliable, high grade cycles (e.g., 
Bradenburg et al. 1998, SB). In this paper, we therefore take 
a fresh look at cycle amplitudes, in light of the new results 
for -Pcyc, using the same high quality cycle sample. 



2. Analysis 

We base our analysis on the Mount Wilson Ca II HK data 
(Bea95). Soon et al. (1994) gives cycle amplitudes for a num- 
ber of inactive stars, and Baliunas et al. (1996) gives them for 
all stars with cycles in Bea95, but the A cyc values for com- 
mon stars in these studies could not be reconciled. We have 
thus resorted, for this initial study, to "eye" estimates of the 
5hk index amplitude AShk (based on plots in Bea95) for 
stars not in Soon et al. (1994). These were converted into a 
Ai?^ K following Noyes et al. (1984; see Table 1). As a test, 
we also made "eye" estimates of AShk for stars in Soon et 
al. (1994) and found (^ C yc(eye)-A cyc (Soon)) = 0.00±0.03, 
indicating good agreement. 

We studied stars chosen by SB for having well defined 
cycles, which in practice meant all cycles ranked "poor" in 
Bea95, and some graded "fair", were ignored. Stars were 
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Table 1. Stellar Data 



HD# 


B-V 


—log 

#HK 




Prot 

[d] 


P c*yc 

[yr] 


Sun 


0.66 


4.901 


0.22 


26.09 


10.0,84 


1835 


0.66 


4.433 


0.15 


7.78 


7.78 


3651 


0.85 


4.991 


0.36 


44.0 


14.6 


4628 


0.88 


4.852 


0.38 


38.5 


8.60 


10476 


0.84 


4.912 


0.38 


35.2 


9.60 


16160 


0.98 


4.958 


0.32 


48.0 


13.2 


18256 


0.43 


4.722 


0.16 


3.0 


6.8 


20630 


0.68 


4.420 


0.14 


9.24 


5.6 


26913 


0.70 


4.391 


0.13 


7.15 


7.6 


26965 


0.82 


4.872 


0.38 


43.0 


10.1 


32147 


1.06 


4.948 


0.42 


48.0 


11.1 


76151 


0.67 


4.659 


0.16 


15. 


2.52 


78366 


0.76 


4.608 


0.27, 0.19 


9.67 


12.2, 5.9 


100180 


0.57 


4.922 


0.17, 0.07 


14. 


12.9, 3.6 


103095 


0.75 


4.896 


0.27 


31.0 


7.30 


114710 


0.57 


4.475 


0.21, 0.12 


12.35 


16.6, 9.6 


115404 


0.93 


4.480 


0.16 


18.46 


12.4 


136202 


0.54 


5.088 


0.21 


16. 


23. 


149661 


0.82 


4.583 


0.35,0.15 


21.07 


16.2, 4.0 


152391 


0.61 


4.448 


0.28 


11.43 


10.9 


154417 


0.57 


4.533 


0.17 


7.78 


7.4 


156026 


1.16 


4.662 


0.37 


21. 


21. 


160346 


0.96 


4.795 


0.44 


36.4 


7.00 


165341 


0.86 


4.548 


0.54, 0.12 


19.9 


15.5, 5.1 


166620 


0.87 


4.955 


0.30 


42.4 


15.8 


187691 


0.55 


5.026 


0.11 


10. 


5.4 


190007 


1.17 


4.692 


0.30 


28.9 


13.7 


190406 


0.60 


4.797 


0.34, 0.15 


13.94 


16.9, 2.6 


201091 


1.18 


4.764 


0.32 


35.37 


7.30 


201092 


1.37 


4.891 


0.21 


37.84 


11.7 


219834B 


0.91 


4.944 


0.29 


43.0 


10.0 



*data for a secondary cycle listed second where appropriate 



Some of the scatter in Fig. 1 is also intrinsic. If stellar 
cycles are like the Sun, activity will be restricted in latitude. 
Stars with different inclinations i will then exhibit different 
apparent A cyc (Radick et al. 1998; Knaack et al. 2001). Anal- 
ysis of the models in Knaack et al. (2001), however, indicates 
that 77% of the measured A cyc should lie within ±15% of the 
mean, with only 23% (those with i < 39°) will range from 
-15% to -46% of (A cyc ). Furthermore, stars more active than 
the Sun are expected to have a wider latitude distribution of 
activity (e.g., Schrijver & Title 2001), and thus should show 
less scatter due to varying i. Thus the effects of i should be 
relatively small, and strongest stars with i? HK ^ -Rhk(©)- 
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Fig. 1. Ai?^ K vs. (i? HK ) f° r trie stellar sample; symbols in- 
dicate F (A), G (o, the Sun's is doubled), and K (box) stars, 
filled symbols are more active stars, size oc P cyc and A cyc 
"quality", dotted lines connect secondary P cyc (x) in some 
oc (i?HK)°' 77 (solid line), implying 

-0.23 

HK/ 



stars. We find Ai? HK 



weighted by cycle grade, with an additional reduction applied 
if Acyc (eye) was unusually difficult to determine. Other stel- 
lar properties were also taken from the list compiled in SB. 
We estimated that the amplitude of the solar secondary P cyc 
(the Gleissberg cycle) at w half of the primary (11 year) A cyc . 



3. Results and Discussion 



We first investigated how Ai? HK depends on (-R HK ) itself 



(Fig. 1). We find AR[ 



HK 



HK/ 



\0.77 



(cr = 0.17 dex; fitting 
the primary -P C yc only). This implies that A cyc oc (-Rhk) "'^ 
- fractional cycle amplitudes decrease somewhat with in- 
creasing activity. 

Some of the scatter about the fit can be explained by an 
additional dependence on B — V color: Ai?^ K is larger in K 
stars than G or F at fixed R' UK . This is seen more clearly by 
plotting A cyc vs. B — V color (Fig. 2), which shows an steady 
increase in A cyc from F through G stars, until a maximum is 
reached by the mid K stars (B — V ~ 1.0). The averages 
by spectral type are (A cyc (F)} = 0.17, (A cyc (G)) = 0.22, 
and (A cyc (K)) = 0.35. This behavior may be the result of a 
dependence of cycle amplitude on fractional convection zone 
depth, up to some limiting value in mid-K stars. 
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Fig. 2. A cyc = Ai? HK /(P^ K ) vs. B - V (symbols as in 
Fig. 1), showing an increasing A cyc with decreasing effective 
temperature in F and G stars, reaching a maximum in mid-K 
stars. 



In contrast to the dependence on B — V, A cyc shows lit- 
tle dependence on rotation, whether expressed as ft, or (see 
Fig. 3) as Po _1 = Tc/ P m t (where tq is the convective 
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turnover timescale; Noyes et al. 1984). Since (R'n K ) shows 
a clear dependence on Ro^ 1 (Noyes et al. 1984), this result 
implies AR^ K has a similar dependence on rotation. Indeed, 
for i?o _1 < 1.5, Ai?y K oc Ro^ 1 (above this i?o _1 , there is 
evidence for cycle amplitude saturation). 

Since A cyc < 1, one might infer that only portion of R' UK 
derives from a cycling dynamo. This suggests that the non- 
cycling (small-scale) component of the dynamo is prominent 
and thus (from Fig. 3) must have a similar dependence on 
rotation as the cycling component. However, this is at odds 
with the lack of any strong rotational dependence for -R^k 
in very inactive, "flat" activity stars where activity is likely 
fueled by a small-scale dynamo alone (Saar 1998). Thus, a 
more likely possibility is that there is significant temporal 
overlap in Phk between cycles, reducing the apparent ampli- 
tude A cyc . In this scenario, the cyclic dynamo dominates the 
rotational dependence of R' HK , despite the apparently small 
A cyc . Enhanced cycle overlap with increasing (Phk) would 
also explain the decrease in A cyc with (Pjjk) 0^8- !)■ 
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Fig. 3. A cyc = AR' HK /(AR' HK ) vs. Po" 1 = r c /P rot (sym- 
bols as in Fig. 1), showing no clear relationship. 



-0.2 



-0.4 



-0.6 



-1.0 



-1.2 



- 


A 

°/' 


1 □□ 

□ D r : 

I7r n 

■/a : 












• - 









-3.2 -3.0 -2.8 -2.6 -2.4 -2.2 -2.0 -l.f 



log 



o/n 



Fig. 4. A cyc = AR' HK / (AR'^} vs. uj cyc /n (sym- 
bols as in Fig. 1), showing an inactive (marked "I"; 

oc (cj cyc /f2) ' 66 ) and an active (marked "A"; A cyc oc 
(w cy c/0) - 85 ) branch. 
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Fig. 5. Same as Fig. 4, but including secondary P cyc - An 
inactive (marked "I"; A cyc oc (ut cyc /il) - 68 ) and an active 
(marked "A"; A cyc oc (uj cyc /Q)°- 8S ) branch are indicated. 



Study of the dependence of A cyc on P cyc with our dataset 
reveals a complex situation. Most inactive stars (defined as 
logP^ K < -4.75), combined with a few active ones, trace 
out a relation similar to that found by Baliunas et al. (1996), 
while most active ones show a similar, nearly parallel rela- 
tion offset at higher A cyc . Specifically, 20 stars (14, or 70% 
of them inactive) show A cyc oc (cj cyc /il) ' 66 (a = 0.082 
dex), and 9 stars (7, or 78% of them active) exhibit A cyc oc 
(w cy c/ri) - 85 (a = 0.081 dex). Only two lower quality de- 
tections cannot be assigned to one of these branches. This 
"branched" structure is similar to that seen between Lo cyc /Vl 
and RO' 1 or R' HK by Brandenburg et al. (1998), SB, and Saar 
& Brandenburg (2001) (and first noted by Saar & Baliunas 
1992 and Soon et al. 1993). There are some differences in 
branch membership, though; for example, active G stars HD 
1835, 20630, and 26913 lie on the active ("A") branch in SB, 
but the inactive ("I") branch here (Figs. 4, 5). 

If we also consider secondary P cy c, we find three lie on 
these branches, while four do not. Including the secondary 
P cyc in the branch fits causes little change: for the "I" branch, 
Ac y c oc {uj cyc /Q)°- 6S (a = 0.082 dex) while for the "A" 



branch, A cyc oc (w C yc/^) ' 88 (a = 0.076 dex). Four of 
the stars show A cyc (primary )/A cyc (secondary) = 2.26±0.18, 
suggesting there may be some "preferred" amplitude ratios. 
All four secondary cycles "unmatched" to a branch lie at low 
A C yc and high uj cyc / and may indicate a third branch - more 
data is needed to confirm this. 

Thus, it appears that A cyc and P C yc are related in stars, 
just as they are in the Sun. Study of stars indicates the rela- 
tionship is multivalued, and depends on rotation. 

The present analysis is limited due to the saturation of 
Ca II emission for Ro^ 1 3, which likely suppresses and 
obscures the visibility of cycles for more active stars. Cycle 
overlap is also a concern. A logical next step would be to in- 
vestigate photometric cycle amplitudes, which in active stars 
are due to spots (rather than the plage/network seen in Ca II). 
Photometric cycle amplitudes saturate at considerably higher 
Ro^ 1 (Messina et al. 2001), permitting study of A cyc in much 
more active, faster rotating stars. Ideally, some normalized 
quantity like the fractional luminosity amplitude AL/L or 
the starspot filling factor fs (e.g., from TiO measurements) 
should be used. 
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